TOPOLOGICAL ANALYSIS OF HIGH RESOLUTION CMB MAPS 
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We report the development of numerical tools for the topological analysis of sub-degree res- 
olution, all-sky maps. Software to be released in the HEALFAST (V0.9) package defines 
neighbour relationships for the HEALPIX tessellation of the sphere. We apply this routine 
to a fast extrema search which scales strictly linearly in the number of pixels, N p . We also 
present a highly efficient algorithm for simulating the gradient vector and curvature tensor 
fields "on-the-fly" with the temperature map, needing only of order N p log 2 N p more opera- 
(T) j tions. 

1 Introduction 

£f~) \ The Microwave Anisotropy Probe (MAP) and the Planck Surveyor satellite missions promise to 

produce high-resolution full sky maps of the Cosmic Microwave Background (CMB) within the 
next decade. These maps will contain a great deal of cosmologically relevant information and 
, therefore present a tremendous opportunity for doing high-precision cosmology. Owing to the 

unprecedented wealth of data, this opportunity brings with it the challenge of developing tools 
Q H i which help extract this information in manageable processing time. 

Much attention has been given to the problem of estimating the angular power-spectrum 
of CMB anisotropies. To illustrate the difficulties one is facing, we note that even choosing the 
fastest presently available algorithm and allowing for parallelisation and technological advance, 
the estimation of the angular power-spectrum C\ of CMB anisotropies will take > 6 years0 for 
future million pixel maps. 

While the C\ spectrum is the statistic of choice for a large class of cosmological models 
where the primordial fluctuations are Gaussian distributed, the possibility remains that another 
physical mechanism produced structure in the universe. In this case the CMB maps will almost 
certainly be non-Gaussian and contain more information than merely the power on different 
angular scales. Failing this, non-linear gravitational effects and foregrounds will produce de- 
tectable non-Gaussian signals. 

Even for Gaussian maps, it has been shown that extrema statistics capture information about 
cosmological parameters and their theory is well developed i'iSi. The importance of studying 
these and other topological properties of CMB temperature and polarisation has been noted by 
several authors! Extrema have the advantage that they dominate the noise and are therefore 
very robust. As a consequence they have been among the first scientific results reported from 
recent CMB observations (one of them even at this conference 0!). 

The aim of this talk is to present tools which allow the implementation of a wide range of 
data analysis and processing techniques on the sphere. 



"Full size, full colour images can be downloaded from http:/ /www. tac. dk/^wandel\ 
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Figure 1: Quadtree pixel numbering. The coarsely pixelised coordinate patch on the left consists of four pixels. 
Two bits suffice to label the pixels. To increase the resolution, every pixel splits into 4 daughter pixels shown on 
the right. These daughters inherit the pixel number of their parent (boxed) and acquire two new bits to give the 
new pixel number. Several such coordinate patches are joined at the boundaries to cover the sphere. All pixels 
carry a prefix (here omitted for clarity) which identifies which patch they belong to. 

2 Tools 

2.1 Neighbour relationships on the sphere 

Finding neighbours of a given pixel is essential for performing many local operations on digitised 
data, such as finding extrema, saddle points and zero crossings, finite differencing, extraction of 
patches, real-space filtering etc. While next-neighbour relationships are trivially defined on a 
lattice in the plane, on the sphere it is necessary to match distinct coordinate patches. Further, 
in a spherical topology the number of next -neighbours is not the same for all pixels. 

We base our methods on the HEALPIX tessellation of the sphered. However, even when a 
pixelisation of the sphere is given, we must still decide how to number the pixels. We choose 
the hierarchical, quadtree numbering scheme illustrated in Figure |]. This scheme is well-known 
and was used for the COBE pixelisation. Its main advantages for finding next-neighbours are 
that 

• a neighbour can be found by direct operation on the binary representation of the pixel 
number which is computationally very cheap (finding all next-neighbours of 12 x 10 6 pixels 
takes less than 10 seconds]]); 

• the algorithm is manifestly resolution independent owing to the self-similarity of the num- 
bering scheme; 

• once an algorithm exists that finds next-neighbours, it can easily be extended to find a 
whole neighbourhood of a pixel by first degrading, then finding next-neighbours and then 
prograding. 

A disadvantage is that in its implementation the boundary matching of coordinate patches 
has to be customised for the chosen pixelisation scheme. But even here the symmetries of the 
quadtree are useful. For example, if two patches like the one on the right of Figure [l] were joined 
side by side, neighbours of pixels at the joining edge can be found by inverting all even bits of 
their pixel number, e.g., 4 = OIOO2 «-> 14 = 11102- This and analogous prescriptions generalise 
to arbitrary resolution. 

As an example application, we show in Figure |2| the theoretical prediction and our numerical 
results for an extrema search in a "standard" CDM and an open model with Vt = 0.2 (Hq = 

6 All timings are "wall clock" for a single SGI RfOOOO 194MHz processor. 
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Figure 2: Probability density for finding a maximum (minimum) above (below) a certain excursion level. The 
minima distributions have been flipped over to coincide with the maxima distributions. Open dots are found 
from one realization of an SI = .2 map, filled dots are for standard CDM. Lines show the theoretically expected 

distributions (solid: SCDM, dashed: OCDM). 



50km/s/Mpc, = 0.05 for both, Ofwhm = 10' and zero noise). While it is clear that these 
probability distributions contain much less information than the anisotropy spectrum (they form 
a 3 parameter family if the total number of pixels is counted), the distributions shown clearly 
distinguish the two models. We stress that this is a realistic situation, since we used data only 
from one map. 

This example do not even begin to explore the possibilities of these methods. By just 
measuring the probability distribution of excursion of an extremum from the mean, we are 
throwing away valuable scale and phase information. A more careful analysis will reveal more 
about the specific advantages of topological analyses of CMB maps. 



2.2 Map simulation and computation of derivative fields 

The derivative fields of a map are important for the study of the properties of its station- 
ary points. A Fast Spherical Harmonics Transform code exists for the synthesis of Gaussian 
HEALPIX temperature maps with a given C\ spectrum!. A priori, generating the derivative 
fields would take just as long for each component as the map generation itself. The time limiting 
factor for map generation is the calculation of the associated Legendre polynomials P ; m (cos#) 
for each ring of constant co-latitude 0. We exploit this fact by calculating the derivative fields 
at the same time as the temperature map and exploiting the following differential relations: 



m \ 



^Y lm {6,4>) = -mcot 9Y lm (0, 0) - + 1) - m(m - l)Y lm ^(9, tfe* (i) 

In this way, only five more Fast Fourier Transforms have to be computed for each ring on the 
sphere, with the total added cost scaling as iVplog 2 N p . For a map with N p = 12.6 x 10 6 pixels, 
this amounts to an added 15 minutes of wall clock computation time for all five derivative 
components compared with 20 minutes for the temperature alone. In Figure ^ we display 
temperature anisotropics, the square norm of the gradient vector |VjT| 2 = si ^ ig (-j^T) 2 + ( J^T) 2 

and the trace of the curvature tensor Tr(Ajj) = ]^?q -§j?T + -§mT, for the standard CDM model. 
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Figure 3: From left to right: the temperature, the square of the gradient vector and the trace of the curvature 
tensor for a standard CDM sky smoothed at 5 degrees for easier visualisation. 



3 Conclusions 

In this talk we presented tools which are generally applicable for the analysis of topological 
properties of maps on the sphere such as the detection of extrema, saddle points and zero 
crossings of a function, even though they can be used more generally for any kind of local 
analysis of spherical data sets. The computational time required for all sky computations scales 
strictly as the number of pixels in the map, N p . Further, we extend existing simulation tools for 
Gaussian fields to generate not only temperature maps but also the associated gradient vector 
and the curvature tensor fields, with the extra computational cost scaling as iVplog 2 N p . 
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